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Abstract 



The exact perturbation approach is used to derive the elementary correlation lengths 
^i and related mass gaps mi of the two-dimensional dilute Al lattice model in 
^ . regimes 1 and 2 for L odd from the Bethe Ansatz solution. In regime 2 the ^3 

model is the E^ lattice realisation of the two-dimensional Ising model in a mag- 
netic field at r = Tc. The calculations for the A3 model in regime 2 start from 
CN| ■ the eight thermodynamically significant string types found in previous numerical 

, studies. These string types are seen to be consistent in the ordered high field limit. 

Q«^ \ The eight masses obtained reduce with the approach to criticality to the E^ masses 

"^ ' predicted by Zamolodchikov, thus providing a further direct lattice determination 

of the £"8 mass spectrum. 
Oh! 



1 Introduction 



There is an intimate relationship between conformal field theory, integrable 
field theory and integrable lattice models in statistical mechanics [1]. In par- 
ticular, massive integrable field theory can be considered as conformal field 
theory perturbed by some scalar relevant operator. The perturbed Hamilto- 
nian is 



H = H, + X f(^{r)d\r). (1.1) 



The canonical example is the Ising model. The $(1,3) perturbation is thermal 
and introduces a single correlation length into the system. The off-critical 
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system is described by scalar combinations of massive Majorana fermions of 
mass proportional to the inverse of the correlaton length. Here the integrable 
field theory can be considered as c = | conformal field theory perturbed by 
the energy density $(1,3) having the conformal dimensions (|, |). 

In a remarkable advance, Zamolodchikov [2,3] considered the $(1,2) magnetic 
perturbation and showed that there are a number of nontrivial local integrals 
of motion and thus an integrable field theory. In this case the c = | conformal 
field theory is perturbed by the spin operator of dimension {jq,jq)- Zamolod- 
chikov then conjectured the S'-matrix and mass spectrum of this field theory. 
The masses coincide with the components of the Perron-Frobenius vector of 
the Cartan matrix of the Lie algebra Eg. They are 

m2/mi=2cosf =1.618 033... 

m3/mi=2cos§ =1.989 043... 

m4/mi = 4 cos f cos If =2.404 867... 

mg/mi =4cosf cosff =2.956 295... (1-2) 

mg/mi = 4 cos f cos § =3.218 340... 

rm/mi = 8 cos^ | cos |f = 3.891 156 .. . 

ms/mi = 8 cos^ | cos ff = 4.783 386 .. . 



Very soon after, numerical tests were performed to check that these masses 
were indeed present in the lattice model, namely the two-dimensional Ising 
model in a magnetic field at T = T^. The first few masses were convincingly 
observed in the related quantum Ising chain in a magnetic field via finite-size 
diagonalisation [4-6] and the truncated fermionic state space method [7]. 

Zamolodchikov also pointed out the possibility of an integrable off-critical lat- 
tice model corresponding to the integrable magnetic perturbation of the field 
theory. In a further development, an integrable lattice realisation of the Eq 
Ising model was provided by the dilute A^ model [8]. In this model the elliptic 
nome plays the role of magnetic field. The calculation of the bulk free energy 
of the dilute A3 model in the appropriate regime gives the magnetic Ising ex- 
ponent (5 = 15 [8], which also follows from the calculation of the local height 
probability [9]. The expected Ising magnetic surface exponent (5^ = — y fol- 
lows from the excess surface free energy [10]. The study of the thermodynamics 
of the dilute A3 model revealed the entire E^ mass spectrum in the scaling 
limit [11]. In particular, the resulting integral equations in the thermodynamic 
Bethe Ansatz calculations are those discussed earlier based on the Lie algebra 
Es [12]. The Eg structure in the dilute A3 model in regime 2 has also been 
established by expressing the one-dimensional configuration sums appearing 



in the local height probability [9] in terms of fermionic sums which explicitly 
involve the Eg root system [13]. This "fermionic sum = bosonic sum" expres- 
sion yields the Eg Roger s-Ramanuj an identity for the Xi i Virasoro character 
[14]. The Fourier transform results for the single particle excitations [11] have 
been inverted to obtain analytic expressions for the excitation energies of the 
eight quasiparticles in the Hamiltonian version of the model [15]. 

Including the groundstate, the thermodynamic Bethe Ansatz relied on the 
input of nine sets of thermodynamically significant string solutions of the 
Bethe equations. These string types were checked numerically as far as possible 
at criticality [1 1] . The stability of these strings types away from criticality has 
been examined in detail [16,17]. One of the string types, associated with the 
mass 7714, was seen to differ from that proposed in [11]. String distributions 
were observed up to mass m^, with a number of elementary excitations up 
to mass rrii + m^. However, root distributions associated with the masses 
rriQ to TTig were not observed, presumably because of the inherent numerical 
difficulties. 

In this paper we explicitly derive the inverse correlation lengths, and thus the 
mass gaps, of the dilute A^, lattice model. Q We use the string solutions given 
in [11,17] and apply Baxter's exact perturbative approach [19], as used for 
example in the calculation of correlation lengths in the cyclic solid-on-solid 
(CSOS) model [20]. In particular, we obtain the mass gaps given in (1.2) as 
criticality is approached. 

The outline of the paper is as follows. The dilute A^ lattice model is defined 
along with the corresponding Bethe equations in Section 2. The bulk free 
energy is derived via the exact perturbation approach in Section 3. The eigen- 
value expression for the leading excitations in regime 1 is derived for L odd in 
Section 4. The eigenvalue expressions in regime 2 for L = 3 associated with 
the eight Eg masses are derived in Section 5. These results are collected to- 
gether in one formula in Section 6, where the corresponding correlation lengths 
and masses are given. The paper concludes with a discussion of the results in 
Section 7. Some intermediary results are given in the Appendices. 



2 The dilute A^, lattice model 



The dilute A^, model is an exactly solvable, restricted solid-on-solid model 
defined on the square lattice. At criticality, the model can be constructed 
[8,21] from the dilute 0{n) model [22,23]. Each site of the lattice can take 
one of L possible (height) values, subject to the restriction that neighbouring 
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sites of the lattice either have the same height, or differ by ±1. The off-critical 
Boltzmann weights of the allowed height configurations of an elementary face 
of the lattice are [81 
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The crossing factors S{a) are defined by 
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and 'i9i('u), '*?4('u) are standard elliptic theta functions of nome p, 



n=l 



^^{u)=di{u,p) = 2p^/^smu n (l- V"cos2m + p^") (l-p2"j, (2.3) 

(2.4) 
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In the above weights the range of the spectral parameter u and the variable 
A are given by < -u < 3A with 



A = - TT, 

r 



(2.5) 



where r = 4(L + 1), with s = L in regime 1 and s = L + 2 in regime 2. We do 
not consider the other regimes here. The thermal Ising point occurs in regime 
1 with L = 2 and the magnetic Ising point occurs in regime 2 with L = 3. 

The row transfer matrix of the dilute A models is defined on a periodic strip 
of width A^ as 
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(2.6) 



where {a} is an admissible path of heights and aN+i = ai, b^+i = &i- For 
convenience we take N even. 

The eigenvalues of the transfer matrix are [11,24,25] 



A(«) 



-UJ 



+ 



^i(2A-m)^i(3A-m) 



-|7V TV 



^i(2A)^i(3A) 



n 



'i9i(-u — Uj + A) 



^i(m)^i(3A-m) 
^i(2A)^i(3A) 



-iTV TV 



n 



f=i^i{u-Uj - A) 

'i9i(-u — Uj) diiu — Uj — 3A) 
\ di{u — Uj — A) di{u — Uj — 2A) 



+ UJ~ 



§i{u)di{\-u) 
^i(2A) ^i(3A) 



N TV 

n 



§i{u-Uj -AX) 
\'di(u — Uj — 2A) 



(2.7) 



where the A^ roots Uj are given by the Bethe equations 
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and tu = exp( i 7r£/(L + 1)) for £ = 1, ... , L. 

There are several methods at hand to calculate the correlation length. Here 
we apply the perturbative approach initiated by Baxter [19,20]. For L odd 
this involves perturbing away from the high magnetic field limit aX p = 1. We 
thus introduce the conjugate variables 
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where nome p = e ^. The relevant conjugate modulus transformations are 



where q = e~^ '^ and 
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In the ordered limit (p — > 1 with ujt fixed) the Boltzmann weights for L odd 
reduce to 
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The function H{d,a,b) is given explicity in [9], being required for the calcu- 
lation of the local height probabilities. In this limit the row transfer matrix 
eigenspectra breaks up into a number of distinct bands labelled by integer 
powers of w. In regime 1 there are |(L + 1) ground states and in regime 2 
there are |(L — 1) ground states, each with eigenvalue Aq = 1. The bands of 
excitations are relevant to the calculation of the correlation lengths. 



The number of states in the w band is ^{L — 1)N in regime 1 and ^{L — 3)N in 
regime 2. These correspond to introducing in all but one of the ground state 
paths {a} a single non-ground state height, in any position. In particular, note 
that there are no excitations in the w band for L = 3 in regime 2. Thus for the 
magnetic Ising model the leading excitations are in the w"^ band. These are 
harder to count, arising from a variety of both single and multiple deviations 
from ground state paths. However, we observe numerically that (apart from 
when N = 2) there are 4N states in this w"^ band. 

In our numerical investigation of the transfer matrix eigenspectrum we asso- 
ciate a given value of i with each eigenvalue by comparing the eigenspectrum 
at criticality {p = 0) with the eigenspectrum of the corresponding 0(?t,) loop 
model for finite N.^ Each eigenvalue can then be tracked to the ordered 
limit. In this way the band of largest eigenvalues is seen to have the values 
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value of i for each ground state. 



Setting Wj = e 27r%/e^ ^]-^g eigenvalues (2.7) can be written 
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The Bethe equations (2.8) are now 
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We are now ready to investigate the ordered limit. 
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^ Strictly speaking we compare with the eigenspectrum of the corresponding vertex 
model with seam oj. 



3 Free energy 

The calculation of the largest eigenvalue in the thermodynamic limit proceeds 
from the a; ^ limit with w fixed in a similar manner to that for the eight- 
vertex [19] and CSOS [20] models. We make repeated use of the properties 

E{z,p) = E{p/z,p) = -zE{z-\p). (3.1) 

Assuming that the roots Wj are on the unit circle, as observed in our numerical 
calculations, only the first term in the eigenvalue expression (2.14) survives in 
this limit, with 

Ao-cc;(«;i...w^)("-2s)/r^ ^3^2) 

The Bethe equation (2.15) gives 

w^ + u-\wi...WNf'''' = Q, (3.3) 



where we write w = Wj. We consider this equation for all complex w and 
equate the left hand side to Y[f=i{w — Wj) so that 

(wi . . . wn) = u;-\wi . . . wjv)''/' ^{wi... wjv)^'"''^/' = u;-\ (3.4) 



Recall that we have taken A^ even. Thus from (3.2) the largest eigenvalue is 
Aq = 1 in the limit x -^ 0. Further, we expect the product of the roots to obey 
(3.4) away from a; = 0. Each of the degenerate ground states has a different 
root distribution {wj} on the unit circle, depending on £. 

To perturb about a; = 0, we define the auxiliary functions 
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where A{z) is known, and Fq{w) and Go{l/w), which depend upon the Wj, are 
the unknowns to be found. Then the Bethe equation (2.15) may be written as 



which is again an A^th order equation with roots Wi, . . .Wat, so that the left 
hand side may be equated to 
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These products can be expressed in terms of the auxiliary functions, 
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where we have made use of (3.4). When the above solutions for J^q{w) and 
Qq{1/w) are used, we can show that all three terms are identical, so that we 
finally obtain 
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This is in agreement with the previous calculations via the inversion relation 
method [8,9]. The singular behaviour follows as 
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where the exponents are given by [8,9] 
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These include the Ising values a = for L = 2 in regime 1 and 5 = 15 for 
L = 3 in regime 2. 



4 Excitations in regime 1 



In regime 1, we observe numerically that the leading eigenvalue in the w band 
has £ = |(L+1) + 1. The corresponding root distribution has A^ — 1 roots 
on the unit circle and a 1-string excitation located exactly at Wn = —x"^. In 
general we assume that the 1-string excitations are located at wn = bx^ with 
\b\ ~ 1 and the remainder on the unit circle, i.e. Wj = aj for j = 1, . . . A^ — 1. 
The Bethe equations split into two sets, the first of which is 
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for k = 1, . . . ,N — 1, where we have defined HTLi Oj = ^m- For k = N the 
other equation is 
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which taken together give a'^ = 1, where a = ak- This equation is of order A^, 
but recall there are only A^ — 1 of the a^'s. The extra degree of this equation 
defines a hole at a^. Further, from (4.3) and (4.4) we also have 

cu-\Aj,_^bf'/' = Ar, = Ajv-iajv = -1. (4.5) 
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The Bethe equation for k = 1, . . . ,N — 1 can now be written 
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We are thus led to define 
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where J-'o and ^o are the functions used in the previous Section. 

We can now construct the eigenvalue, which in terms of the auxiliary functions 
reads 



A{x^'-'^'w)A{x'^'-^'w)A{x'^'/w)A{x^'/w) 



X 



A{x'^')A{x^')A{x^'-^')A{x^'-^') 

R(^X^r^^'w)S{l/x^'-^'w) 



+ 



R{x^'w)S{l/x^'w) 

X{x''^^'w)Y{l/x'"'^'w) ^ ' ^ ' ' 
A{w)A{x^''-^'w)A{x'^'' lw)A{x^' jw) 
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+ 



R{x'^'-*'w)R{w/x'^')S{l/x'^''-^'w)S{x'^'/w) 
^ R{w)R{x^^-^'w) S{l/w)S{l/x^''-^'w) 

X{x'-^'w)X{w/x'') Y(l/x'-^'w)Y{x''/w) 
^ X{x'-*'w)X{w/x^+^') Y{l/x'-^'w)Y{x'+'^'/w) 

J-'{w/x'^'^)Q{x'^^/w) 
x"^'-" A{w)A{x'^''/w)A{w/x'^')A{x'^''+'^'/w) 
AN-ih A{x^')A{x^')A{x'^''~*')A{x'^^-^') 
R{w/x^')S{x^'/w) 
R{x^'-^'w)S{l/x^'-^'w) 
X{x'^-^'w)Y{l/x''~^'w) 1 

'X{w/x''+^')Y{x''+^'/w) T{w/x*')g{x*'/w) ' 



X 



(4.17) 



Making use of the solutions (4.15) and (4.16), the identities in Appendix A, 
and (4.5) on the prefactors, we can show that this simphfies to 



A/Ao = wAn-1 



E{x'^^w/a]^^ x^'^^)E{x'^^w/a]\f, x^'^^) 



(4.18) 



If both the hole a^, and b are equal to —1, then we have 



A 



12s > 



W 



E{—x'^-yw, a;^^*) E{—x'^^/w, x 
E{-x^' w, x^"^') Ei-x^"" w, x^^') 



(4.19) 



At the isotropic point w 



x^^ (u 



3A/2) this reduces to 



A 

a; 



X 



E\-3fy^ 



^,(iL,p-/6A) 



■12 
'Stt 



^4(ff,r/^^), 



(4.20) 



We can also express the second Bethe equation (4.2) in terms of the auxiliary 
functions. This yields an equation for the parameters b and qn of the form 



A{x'-^'b)A{x''+^yb) _ 

A{x'-^'/b)A{x'-+^'b) ~ 

R{x'-H)R{x''+^'b)S{l/x''-H)S{l/x''+H) 

R{x'-^'b)R{x'-+^'b)S{l/x'-'^'b)S{l/x'-+^'b) 

F{x''-^'b)F{x'+^'b)G{l/x''~^'b)G{l/x''+^'b) 

^ F{x^~^'b)F{x^+^'b) G{l/x^-^'b)G{l/x^+^'b) ' 



(4.21) 
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Since this cannot be expressed in terms of T and Q ^ we have to re-construct 

OO CO 

F(a) = n T{x^'''-^'^'d), G{\la) = [] Qil/x^^'-^'^'a), (4.22) 



e=o i=o 



from the definitions (4.13). Making use of identities hsted in Appendix A we 
finally obtain a result, 



r-4:s} 



E{x^h/ai^, x'^ 



in which neither the original nome x^^ nor x^^^ appears. This is obviously 
satisfied by ajy = b. 



5 Excitations in regime 2 for L = 3 



For L = 3 in regime 2 extensive numerical investigations of the Bethe equa- 
tions (2.8) have led to a conjecture for the thermodynamically significant 
strings [11,16,17]. For example, the leading excitation Ai, corresponding to 
mass mi, is a 2-string with i = 2. However, this state is originally a 1-string 
for small p and small A^. Such behaviour has been discussed in [17]. By track- 
ing this state numerically with increasing p we see that this 2-string is exactly 
located at —x"^^^ in the limit p = 1. Moreover, it lies in the w"^ band of ex- 
citations. The string structure for the other eigenvalues (Aj) or masses {ttlj) 
used in the thermodynamic Bethe Ansatz approach were checked numerically 
as far as possible at criticality [11]. The stability of these strings types away 
from criticality has been examined in detail [16,17]. The conjectured string 
positions are listed in Table 1. With the exception of mass m^, these are the 
string positions given in [11]. Rather than the 5-string, ±8, ±12, 16, proposed 
in [11] we consider the 7-string found in [17]. In the latter study, string distri- 
butions for elementary excitations up to and including mass m^ were located, 
including a number of composite excitations up to mass mi + m^. However, 
the root distributions associated with the masses uiq to mg were not observed, 
presumably because of the inherent numerical difficulties. 

Here we show that each of the tabulated sets of strings provides a consistent 
solution in the ordered and large A^ limits. However, this is also true for the 
proposed 5-string for mass 7724. According to our calculations, both the 5-string 
and the 7-string distributions yield the same mass m^. We shall return to this 
point in the discussion in Section 7. 

The calculations in this Section are necessarily complicated. The reader not 
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Table 1 

String positions and corresponding eigenvalue bands for the eight elementary masses 

rrij of the dilute ^3 model in regime 2. The strings are in units of 7ri/32. 



i 


string positions 


band 


1 


±11 


2 


2 


±5, ±15 


2 


3 


±10, ±12 


3 


4 


±1, ±7, ±13, 16 


3 


5 


±9, ±11, ±13 


4 


6 


±6, ±10, ±14, 16 


4 


7 


±8, ±10, ±12, ±14 


5 


8 


±7, ±9, ±11, ±13, ±15 


6 



interested in the specific details may wish to skip to Sections 6 and 7, where 
the results obtained are summarised and discussed. 



5. 1 Mass iTLi 



We begin the perturbation argument with the structure Wj = aj for j = 
1, . . . ,N — 2 with Wn-i = bix~^^ and Wn = b2X^^. However, from the Bethe 
equations for k = N — 1 and k = N (see Appendix B) we can show that 
61 = 62 = b. The Bethe equation for the other roots a^ = a is then 



—cu 



^ E(x^^a) 



N 



/y E{x%/a)Eix''b/a) ^^' 



E{x^^a/aj)E{x^^aj/a) 



(5.1; 



b^ E{x^a/b)E{x^^a/b) f^^ E{x^Oaj/a)E{x^Oa/aj)' 
where as before, IVjLi % — ^m- In the a; ^ limit this gives the equation 

-' + -{AN-2b'fVb' = 0, (5.2 



,Af- 



(jj 



which is an (A^ — 2)th order equation for A^ — 2 zeros, so that, unlike in regime 
1, there are no holes in this case. Equating this as usual with Ylfs^{a — a 



1)1 



16 



we obtain 



00 



(5.3) 



which we substitute into the other Bethe equations (B.l) and (B.2) in this 
hmit, to give 



UJ 



iAN^2b'r' 



-,2 



{A^-2b'r 

^2N 



6' 



27V 



(5.4) 



Note that the value 5 = — 1, determined numerically for the leading 2-string 
excitation, satisfies this equation. In principle these equations can be solved 
consistently to determine the total number of 2-string excitations. 

To proceed with the calculation we use the auxiliary functions A{w), X{w), 
Y{l/w) defined in (3.5), (4.6) and (4.7) with the appropriate choices s = 
5, r = 16, and define in analogy the "unknown functions". 



7V-2 oo 

j=l fc=0 

N-2 oo 

Gi{i/w)=ii n(i-^'%M, 
j=i k=i 

along with the functions we will actually solve for. 



(5.5) 



gi{l/w) = Giil/w)/Giil/x''w). 



(5.6) 



In fact, although we must define Fi{w) and Gi{w) for each rrii, we can define 



J^,{w) = F,{w)/Fi{x^^w), 
g,{l/w) = G,il/w)/G,il/x'^w) 



(5.7) 



for i 7^ 4, 6. Rearranging the Bethe equation (5.1) and equating the dominant 
terms, as in regime 1, gives 



J^i(a) 
Giil/a) 



Ajx^^a) X{x^^a)X{x^^a) Tijx^^a) 
A{x^^a) X{x^a)X{x^^a) J^i^x^^a)' 

A{x^'^/a) Y{l/x'^a)Y{l/xa) gi{x^^/a) 
A(a;30/a) Y {x"^ / a)Y {x^^ / a) ^i(x20/o) 



(5. 
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m=0 




m=l 



We solve these recursively, using identities for X{w) and Y{l/w) given in 
Appendix A, so obtaining 

^ , , ^ , ,X(x23a)X(x33a) 
^la =-^oH y n ir 53 

1 - x60'"+3ia/6)(l - x60™+39a/6) 
(1 - x60™+9a/6)(l - x60'»+i9a/6) 

^ (1 - a;60™+5ia/5)(i _ a;60™+6ia/6) J ' ^'^ 

and similarly 

1 - x6°'^-3^6/a)(l - x'^^'^-^ift/a) 
(1 - x60'"-96/a)(l - x60'^-i6/a) 

^ (1 _ x60™+96/a)(l - x60™+4i6/a) J ' ^^"^^^ 

We are now in a position to evaluate the eigenvalue (2.15). Expressed in terms 
of the auxiliary functions and using (5.3) it reads 



_ w"^ A{x'^w)A{x^'^w)A{x^^ /w)A{x^^ /w) 
^^ 'W A(x2)A(a;i2)A(x20)A(x30) 

X{xw)X{x^^w) Y{l/x^^w)Y{l/xw) '^ ^^^HV^: w) 
A{w)A{x^w)A{x^°/w)A{x^yw)X{w/x^^)X{w/x^)X{x^^w)X{x^^w) 
^ A(x2)A(xi2)A(x20)A(x30) X(w/x2i)X(x«;)2X(x23w;) 

y(i/x^3w)y(i/x^^w)y(xVw)y(3:"/w) j^i(w)gi(i/w) 

^ F(l/x23«;)F(l/x«;)2r(x2Vu;) J^i(w;/xio)^i(xi7w) 

A{w/x^^)A{w)A{x^yw)A{x^yw) X{w/x^^)X{x^w) 
^ A(x2)A(xi2)A(x20)A(x30) X(w/x3i)X(w/a;9) 
Y{l/x^w)Y{x^yw) 1 



r(x7w)F(x3Vw) J^i(w/a;20)^i(s20/u;) 



(5.11) 



Using (5.9) and (5.10) we find that this result may be written in terms of 
elliptic functions asQ 



^ We point out here that all three terms in the eigenvalue expression are identical, 
once J-'i and Qi are substituted. This feature is common to all masses rrii {i = 
1, ... ,8) when the appropriate functions J^i and Qi have been found and used. So 



Ai _ w^ E{xb/w,x^^) E{x^^b/w,x^^) 
Tq~1^ E{xw/b, a;60) E{x^^w/b, x^o) 

E{x^^w/b, x^o) E{x^^w/b, x^^) 
^ E{x^^/w, x60) E{x^^b/w, x^o) ■ 

Setting 6 = — 1, the leading excitation for L = 
lies in the w"^ band, is thus 



(5.12) 



3 in regime 2, which clearly 



Ai 
Ao 



■w 



E{-x/w, x^^) E{-x^^/w, x^°) 
E{-x w, x60) E{-x^^ w, x^o) 

^(-x^^ w, x^°) E{-x^^ w, x6°) 
^ E(-x3Vu;, x60) ^(-x^Vu;, x^o) 



(5.13) 



At the isotropic point w 



X 



15 



this reduces to 



Ao 



14 ^60\ 



:X 



28 



^2^_a;4^^60)^2^_^14^^. 



Mt„P 



2,16 2,60) ^2 ( 



_^26 ^,60^ 



15' 



'^')^4(i,p^/^^) 



_^4(f,P«/^5)^4(lf,P«/l^). 

where we have written the excitation in terms of the original nome. 



(5.14) 



5.2 Mass 7712 



We begin the perturbation argument with the structure Wj = aj for j = 
1, . . . , N—4: with WAr_3 = bix~^, wn-2 = b2X^, WAr_i = b^x'^^ and w^- = &4X^^. 
From the Bethe equations for k = N — 3, . . . ,N (see Appendix B.2) we can 
show that 6i = 62 = &3 = &4 = b. The Bethe equation for the other roots 
Cfe = a is then 



-UJ 



'' Eix^^a) 



N 



. . ,4,5/8 0^ E{x%/a)E{x%/a) ^^ E{x^'^a/aj)E{x'^'^aj/a) 



n 



(5.15) 



b^ E{x^a/b)E{x^a/b) fj^ E{x^^aj/a)E{x'^^a/aj)' 
In the X — > limit this gives the equation 

a^-4 + l(^^_^54)5/8/^4 _ Q^ (5^^g) 



00 



in discussion of the subsequent masses, we suppress the second and third terms of 
the eigenvalue expression. 



19 



Equating this as usual with n,=i (O' ~ Oj)) we obtain 



u 



-{A^-.by/' = A^.,b\ 



(5.17) 



(which we later apply to prefactors in A2). From the other Bethe equations 
(B.6)-(B.9) in this limit, 



UJ 



-{AM-.tr/A' 



{AN-4by 

^2N 



l2N 



(5.1^ 



We define functions of the roots we wish to find as 



7V-4 00 



j=l fc=0 

N-A 00 



G2{l/w)=\{\{{l 



X 



32fc 



aj/w). 



j=l k=l 



(5.19) 



Treating the Bethe equation (5.15) as before gives, in terms of the functions 
defined in (5.7), the recurrences 



^2(0) 



Q2{l/a) 



A{x^''a) X(x^^a)X(a:^9Q) J^2(x^°a) 
A(a;22a) X{x^a)X{x^a) J^2(a;20a)' 

^(x^Va) Y{l/x''a)Y{l/x''a) g2{x^Va) 



A(x37a) Y{x^^la)Y{x^Tla) ^2(^270) ' 
Solving these we obtain 



(5.20) 



^2{a)=^o{a] 

00 

m=0 



XV9a)X(A2 

X{x^a)X{x^^a) 

J (1 - a;60™+27a/5)(i _ 3.60^+33^/^) 

(1 - a;60'n+3a/5)(i _ a,60m+i3a/5) 

(1 - a;60'»+37a/^)(i _ 3.60^+43^/^) 

(1 - a;60'n+57a/5)(i _ x^^"'+^''a/h) 



X 



and similarly 



(5.21) 



^2(l/a) = ^o(l/a 



r(l/x^a)r(xVa) 
F(a;iV«)>'(a;'V«) 



20 



oo 
m=l 



X 



60m-37 



h/a){l-x' 



60m-27 



h/a) 



'I 



X 



60m- 13 



X 



60m+17 



h/a){l- 
h/a){l 



X 



60m— 7 



h/a) 



X 



60m+27 



h/a) 



(5.22) 



(1 - a;60™-3^/a)(i _ x^'^'^+^/a) J 

We now substitute these into the eigenvalue expression, the first term of which 
gives 



A, 



w^ A{x^w)A{x^^w)A{x^^/w)A{x^^/w) X{x^^w)X{x^^w 



62 A(x2)A(a;i2)A(a;20)A(a;30) 

Y{l/x'^^w)Y{l/x'^^w) 



X 



Y{l/x^^w)Y{l/x'^w) 



X(x%)X(xi%) 

:F2{x^''w)g2{i/x^''w). 



(5.23) 



This in turn gives an expression in elhptic functions analogous to (5.12), but 
we immediately set h = —1, yielding the leading 4-string excitation in the w'^ 
band for L = 3 to be 



A2 

Ao 



■w 



Ei-x^w, x^^) E{-x^yw, a;6°) 
E{-x'^ w, a;60) E{-x^'^ w, x^^) 

E{-x^'^ w, x^^) E{-x^'^ w, x^°) 
^ Ei-x^yw, a;60) E{-x^yw, x^^) 



At the isotropic point w = x this reduces to 



(5.24) 



Ao 



■x 



-10 



^2(_^2^^60)^2(_^8^^60) 



E2(-a;22,a;60)E2( 



.^28 ^60^ 



^4(i.//^^)^4(ff,P^/^^) 
^,(7|^p8/15)^^(nzE,p8/15) 



(5.25) 



5.3 Mass m^ 



We begin the perturbation argument with the structure Wj = aj for j = 
1, . . . ,N — 4 with wn-3 = hix~^^, wn-2 = &2a:^°, wn-i = hsx^^'^ and wn = 
640;^°. From the Bethe equations for k = N — 3, . . . ,N (see Appendix B.3) we 
can show that 61 = 62 = ^3 = ^4 = h. The Bethe equation for the other roots 
is 



—00 



Ejx^ya) 
'' E(x^^a) 



N 



(A^-46^)^/^ 



64 



21 



^ E{x^a/b)E^{x^^a/b)E{x^^a/b) f}^ E{x^Oaj/a)E{x^^a/aj)' 
In the X ^ limit this gives the equation 



(5.26) 



(5.27) 



which we equate the Ihs as usual with n,=i (o ^ o,j) to obtain 



u 



-{A^^.by/' = A^.,b' 



(5.2J 



Using this with the other Bethe equations (B.14)-(B.17) in the x 
we obtain 



CO 



iAN-4by 



l3N 



1. 



limit 



(5.29) 



Note that since A^ is even, b = —1 satisfies this equation. 
We define the appropriate functions for the unknown roots, 



Af-4 oo 

j=l k=0 

N~4 oo 

G3(i/^)=n n(i-^''s/^)' 
j=i k=i 

then from (5.26) come the recurrences. 



(5.30) 



Ajx^'^a) X{x^'^a)X\x^^a)X{x^^a) J^^jx^'^a) 
A(a;22a) X{x'ia)X^{x^'^a)X{x^'^a) T^ix^'^a)' 

A{x^^/a) Y{l/x^a)Y\l/x^a)Y{l/a) ^3(^:^7^) 
A(a;30/a) Y{x^ /a)Y^{x^'^ /a)Y{x^^ /a) g^{x'^'^/a) ' 

The solutions are 



^s(a) 



Gz{l/a) 



(5.31) 



JFofa) =J^n(a 



01 

oo 

m=0 



X(a;^^a)X(x3'a)X(a:34a) 
'x(a;i0a)X(a;i2a)X(x20a) 



a; 



60m+22 



a/6)(l-a;60'"+30a/6)(l-x' 



60m+32 



a/6) 



x60™+8a/6)(l - x60'»+iOa/6)(l - x' 



60m+18 



a/b) 



22 



(1 - a;60™+38a/6) (1 - x^o^^+^Oa/fc) (1 - x'^^^'^+^^a/h) 1 
^ (1 - a;60'n+52a/5)(i _ a;60™+60a/5)(i _ 3,60^+62^/5) j ' (^-32) 



F(l/a)F(l/x2a)F(xVa 




m=l 



r(a;i2/a)F(x20/a)F(a;22/a) 
1 - a;60"^-326/a)(l - a;6°'^-306/a)(l - x'^'^'^-^^ / a) 
(1 - a;60™-i86/a)(l - x^'^'^-^%/a){l - x^^'^-^/a) 



(1 - x'^''"'+^%/a){l - 3;60"^+20&/a)(l - x'^''"'+^% / a) \ 
^ (1 _ a;60m-25/Q) (1 _ a;60m5/^) (^ _ x'^'>"'+%/a) J ■ ^ ■ ^ 

In terms of these functions the eigenvalue may be represented as 

w^ A{x^w)A{x^^w)A{x^yw)A{x^yw) X{x^^w)X{x^^w)X{x^^w) 



A 



3 ■ 



63 A(x2)A(a;i2)A(a;20)A(x30) X(a;2w)X(xi0w)X(a;i2u;) 

Thus, apphcation of the perturbation argument has yielded the leading exci- 
tation in the w^ band to be 



A3_ 3 E{-xyw, a;^°) E{-x^yw, x^^) E{-x^^/w, 3;^°) 
A^~^ E{-x^ w, x60) E{-x^o w, ^60) E(-a;i2 w, x^^) 
E{-x^^ w, x^°) E(-a;^° w, x^^)E{-x'^^ w, x^°) 
^ E(-a;32/u,, a;60) ^(-x^o/u,, a;60)E(-x42/«;, a;60) ' ^ ' '^ 

where we have put b = —1. At the isotropic point w = x^^ this reduces to 



A3 
An 






2 



(5.36) 



5.4 Mass 7714 



We begin the perturbation argument with the 7-string Wj = aj for j = 
1, . . . ,N — 7 with wn^q = bix~^, W7V-5 = b2X, wn^^ = b^x"'^, w^^s = b^x'^ , 
W7V-2 = b^x'^^, tfAf-i = bgx^^ and w^ = bjx^^. From the Bethe equations 
for k = N ~ 6, . . . ,N, we show in Appendix B.4 that 61 = 64 = 65 = a and 
b2 = b^ = be = j3. We further let 67 = b. This feature is different from that 
seen for the previous masses (or will see for ms, m-j, 7n^). However, in the final 
analysis, our answer does not depend upon a or [3. 
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The other Bethe roots satisfy 



—U) 



'' E(x^^a) 



N 



{A^^ja'P^b) 



3,q3n5/8. 



a2/5252 



E{x^a/a)E{x^^a/a)E{x^^a/a) E{x(3/a)E{x^^(3/a)E{x^'^(3/a 



E{xa/a)E{x^^a/a)E{x^^a/a) Eix'^a/ /3)E{x^'^a/ (3)E{x^'^a/(3) 



E{x^b/a)E{x%/a) ^^ E{x^^a/aj)E{x^^aj/a) 
^ E{x^a/h)E{x^a/h) ^ E{x^^aj/a)E{x'^^a/aj) 



In the X — s> hmit this gives the equation 

a^-6 + l(Ajv-7«'/?'&)'/V«'/3'&' = 0, 



UJ 



(5.38) 



which is an (A^ — 6)th order equation for N — 7 zeros, so that there will be a 
hole, aAr_6, in this case, as there was in regime 1. Equating the Ihs of (5.38) 
with YlfSi^{a — ttj) we obtain 



UJ 



(5.39) 



From the other Bethe equations (B.22)-(B.28) in the a; — > limit, 

i(^.W/^»,)^/f^ '-^--')W^y . ,5.40) 



Now (5.38) must be satisfied by each of the a^, j = 1, . . . , iV — 7, 

-AN_^aN-(i, (5.41) 



<-« 



so that forming the product of these (A^ — 7) equations, 

(A^_r)^-^ = -{AM-iaN-.f-' => {a^-ef = -A^^riaN-eY ■ (5.42) 



Taking (5.39), (5.40) and (5.42) together gives 



{b'^ttN- 



\N 



(5.43) 



which is clearly satisfied by 6 = — 1 and Qn-q 
analogy with the previous auxiliary functions. 



-1. We need to define, in 



24 



Xa{w) = n(l - x^'^'^w/a), Ya{l/w) = n(l - a:'''«/u'), (5.44) 

A:=0 fe=l 



oo oo 



Xp{w) = l[{l-x'''w/P), Yp{l/w) = l[{l~x'''p/w), (5.45) 

fc=0 k=l 



oo 



i?4(w) = n (1 - a^^^'^w^/aw-e), (5.46) 



k=0 

oo 



54(l/ti;) = n (1 - x'^'a^-a/w), (5.47) 

fc=i 

and the functions 



7V-6 oo 

j=l fc=0 

N~6 oo 

6-4(1/^)= n n(i-^^^''«i/^)' (5-48) 

j=i k=i 

for the as-yet-unknown roots. In place of the usual definition, we have 

Fii{w) Xa{xw)Xp{x'^w) 



Ta{w) 



F4{X^^W) Xa{x^w)X(^{x^^w) 



G^jl/w) Y^il/xw)Y^il/x'w) 
^^ ' G'4(l/a;i2u;)F,(l/x5u;)F^(l/a;i%)' ^ ' ' 

Rearranging the Bethe equations (5.37) gives recurrences as usual, with 
Aix^'^a) X(x26a)X(a;28o)i?4(a:20a)i?4(x22a) T^ix^'^a) 



TA\a) 



A(a;22a) X{x^a)X[x^a)R^{x^'^a)R4{x^'^a) T^ix^^a) 



^(x^Va) Y{\lx^a)Y{\lx^a)Sl\lx'a)S4{\la) Q^jx^'^la) 
^'^ '""^ Aix^^ja) Y[x^^la)Y[x^^la)S4[x^la)S4[x^^la) Q^ix^'^ja) '^ ' 

which we solve using identities for the various auxiliary functions listed in 
Appendix A to obtain 



^,, ^, ,X(x^'^a)X(x^'^a)RJx^^a 



X{x^a)X{x^'^a) R^ix^^a) 

00 ( (I _ „60m+26aw-| _ 60m+34aw-| _ 60m+36a^ 
X 



n 



n — a;60m+4a\Q _ j.60m+Ua\f^ _ ^e0m+5e a.\ 



'I _ r^60m+A4a\ n _ j,60m+30 a \n _ ^6 



y. ^^ bj__ "iV-6^^ aN-6' I /p- c--|\ 

^l _ 2,60m+66a^ Q _ j,60m+10 a \n _ ^60m+60 a \ I ' ^ > 

\ b' y aff-a '^ aN-6' 
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and similarly 

c n in\-c n / X ni/^!^m^>i jMIZ^ 

X n ' ^ "^^ 



-)(1- 



X 



m=l 



-X' 



.60m-14 b 



-)(1 



■ X 



SOm-eb 



X 



X 



,60m+26 6 



^)(1-X 



,60m-30 CAT-e 



-)(1- 

^)(1-X 
a ' ^ 



SOm-4 b 



X 



60m+20aiV-6 



(I _ x^^™-^*-)(l — a;60ni-10M^\Q _ ^, 



,6nm °JV-6 




The eigenvalue expression is 



(5.52) 



A. 



X 



w^ A{x^w)A{x^^w)A{x'^^/w)A{x^^/w) 
aN-6b^ A{x^)A{x^^)A{x^^)A{x^^) 
X{x^^w) Y{l/x'^^w) 



X{x^w) Y{l/x^w) 
which, in elliptic functions is 



J^4{x^^w)g4{l/x''>w) 



(5.53) 



A4 



w^ E{x%/w, x^^) E{x^^aN-6/w, x^^) E{x^%/w, x^^) 



Ao aN-eb^ E{x^w/h, x^^) E{x^0w/aN-6, x^o) E{x^^w/b, x^o) 
E{x^^w/b, a;60) E{x^^w/aN-6, x^^)E{x^^w/b, x^^) 



X- 



E{x^%/w, a;60) Eix^^aN-a/w, x^^)E{x^%/w, x^^) 
Finally, with h = aAr_6 = — 1, 



(5.54) 



A 



4 _ 3 E{-x^/w, x^°) E{-x^°/w, x^°) E{-x^^/w, x^^) 

To'^ E{-X^ W, X60) E{-X^O W, X60) E(-Xl4 w, X60) 



X 



E{-x^^ w, x^^) E{-x^^ w, x^^)E{-x 



„44 



w,x 



60\ 



E{-x^yw, x60) E{-x*yw, x60)E(-x44/w, X' 



60 ^ 



At the isotropic point w 



X 



15 



this becomes 



A4 
Ao 



^a(^,V^/^^)^a 



20 



■ 12' 



P 



,8/15) ^^ 



^60 



,p' 



,8/15 ^ 



^^(ZE,p8/15)^^(^^p8/15)^^(|ZL,p8/15) 



(5.55) 



(5.56) 



5.5 Mass m^ 



We begin the perturbation argument with Wj = aj for j = 1, . . . , A^ — 6 and 
wn-5 = bix-^^, wn-4 = b2X^^, wn-3 = b^x'^^ , wn-2 = biX^^, wn-1 = b^x-^, 
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wn = b^x^. In Appendix B.5 we show that the bi are equal, and we call them 
b. The Bethe equation for the other roots is 



—u 



E{x^ya) 
'' E(x^^a) 



N 



l,a^ E{x'^b/a)E\x%/a) 
^ ^-6 ) 1,6 E{xTa/b)E^{x^a/b) 



X 



E\x^^b/a)E{x^%/a) ^ E{x^^a/aj)E{x^^aj/a 



n 



E^{x^^a/b)E{x^^a/b) f^^ E{x^Oaj/a)E{x^^a/aj)' 
In the X ^ limit this gives the equation 



UJ 



(5.57) 



(5.5J 



which is an (A^ — 6)th order equation for A^ — 6 zeros, so that no hole needs 
to be considered. Equating the Ihs with 11^1^ (a ^ <^j) gives 



HA^^^br^' = Aj,.,b'. 

CO 



(5.59) 



From the other Bethe equations (B.31) 



u 



iA^-^br^' 



(A 



N-6'^ 



l6\6 



547V 



lAN 



(5.60) 



For the unknown roots we define the functions 



7V-6 oo 

j=l fc=0 

N-6 oo 

G,{l/w)=ll ll{l-x'''a,/w), 
j=i k=i 

which, after rearranging (5.57), obey the recurrences 



(5.61) 



J^'^(a) 



Ajx^^'a) X{x^'^a)X^{x^^a)X^{x^''a)X{x^'>a) J^5(x^°a) 
A{x^^a) X(x7a)X2(x9a)X2(a:iia)X(xi3a) J^^ix^^^a)' 



Ajx^ya) Y{l/x^a)Y\l/x^a)Y^{l/xa)Y{x/a) G,{x^Va) 
^^^ ''^' A{x^'^la)Y{x''la)Y\x'^la)Y'^{x^^la)Y{x^^la)Q^{x'^'^lay ' 

The solutions are 
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T^ia] =Tn{a 



' X{x^a)X{x^^a)X{x^^a)X{x^'^a) 




X 



60m+23 



a/h){l 



X 



60m+29 



a/6)(l 



X 



60m+31 



a/b) 



m=0 



a,60m+7a/5)(i _ a;60'n+9a/5)(i _ 2,60™+!!^/^) 



(1 - x60™+33a/6)(l - x*^o™+37a/6)(l - x^o^+^^a/fe) 

(1 - x60'n+17a/6)(l - x60™+53a/6)(l - a;60m+59a/5) 



;i - x60™+6ia/6)(l - x60™+63a/^) 



(5.63) 



g,{l/a) = go{l/a) 



Y{l/x^a)Y{l/xa)Y{x/a)Y{xya) 
Y{x^ya)Y{x^ya)Y{x^ya)Y{x^ya) 




X 



60m-33 



b/a){l - x^^"''-^^b/a){l - X 



60m-29 



b/a) 



(1 - a;60™^i75/a)(l - a;60™-ii6/a)(l - x^'^"''%/a) 



60m-23 



b/a){l 



X 



60m+13 



b/a){l 



X 



60m+19 



b/a) 



1 — X 



60m— 7 



1 



X 



60m+21 



b/a){l 
b/a){l- 



X 



60m— 3 



b/a){l 



X 



60m— 1 



b/a) 



X 



60m+23 




(1 - a;60™+i5/a)(i _ ^eo'^+^fe/a) 
which we next substitute into the eigenvalue expression 



(5.64) 



A.. 



w'^ A{x'^w)A{x^'^w)A{x^^/w)A{x^^/w) 
V A(a;2)A(a;i2)A(a;20)A(a;30) 
X{x^'^w)X{x'^^w)X{x'^^w)X{x'^^w) 



X 



X{x^w)X{x^w)X{x^^w)X{x^^w) 
Y{l/x'^'^w)Y{l/x'^^w)Y{l/ x'^^w)Y{l/ x^'^w) 
Y{l/x^w)Y{l/x^w)Y{l/x^^w)Y{l/x^^-w) 
x^5(x'°«^)^5(l/a:'°ti^), 



X 



to obtain (with b 



(5.65) 



As 
Ao 



■w 



X 



^E{-x^ /w)E{-x^ /w)E{-x^^ /w)E{-x^^ /w) 
E{-x^w)E{-x^w)E{-x^^w)E{-x^^w) 
E{—x^^w)E{—x^'^w)E{—x*^w)E{—x*^w) 
E{-x^yw)E{-x^yw)E{-x^yw)E{-x^yw)' 



The eUiptic functions are of no me x . At the isotropic point w 
leading excitation reduces to 



(5.66) 



X 



15 



the 



As 
Ao 



^,(§E,p8/15)^^(^^p8/15)^^(l|^^p8/15)^^(Z|,p8/15) 



30 



15' 



(5.67) 
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5. 6 Mass rriQ 



We begin with Wj = aj for j = 1, . . . , A^ — 7 and wj^_q = bix~ , wn~5 = b2X , 
W7V-4 = b^x^^^, W]\f^3 = b^x^^, W]\f^2 = b^x'^, tfAf-i = bgx^, wn = bjx^^. In 
Appendix B.6 we show that bi = b2 = b^ = bQ = bj = b and that 63 = 64 = a. 
This feature is similar to that for m^, again, in the final analysis, our answer 
does not depend upon a. The Bethe equation for the other roots is 



-co 



E{x^yaj 



N 



(A.-,tV)^'^^^ 



X 



E{x^Oa) 

E{x%/ a)E{x%/ a)E{x%/ a)E{x%/ a)E{x^%/ a) 
E{x^a/b)E{x^a/b)E{x^a/b)E{x»a/b)E{x^'^a/b) 



E{x^^a/a)E{x^^a/a)E{x^^a/a) ^^ E{x^^a/aj)E{x^^aj/a) 
^ E{x^Oa/a)E{x^^a/a)E{x^Oa/a) f}, E{x^Oaj/a)E{x^^a/aj)' 

In the X ^ limit this gives 





N~7 



-(AN-jb^a^f/yb'^a 

UJ 



(5.6J 



(5.69) 



as the equation for the (A^ — 7) roots. Equating this with 117=1 (ct — dj) gives 



-{AN-jb^'a'f'^ = A^^rb^^a. 

UJ 



(5.70) 



From the other Bethe equations (B.36) we obtain 



i{Aj,_rb'a'M' 

UJ 






LAN 



1. 



(5.71) 



We define 



7V-7 00 

j=l k=0 

N-7 00 

i=i fc=i 
and in place of the usual definition. 



(5.72) 



^r(w) 



Fe(x^^w) XJx^w) ' 
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where X^iw) and Fall/"^) have the same definition (5.44) as they did for rui. 
We solve the recurrences 



A(a;iOa) X(a;i8a)X(x24a)X(a;26a)X(a;28a)X(x30a) J^6(a:'°a) 
'^^"^ ~ A(a;22a) X(a;2a)X(a;4a)X(a;6a)X(a;8a)X(a;i4a) JS(a;20a)^ 



i\yx — ja) I yx"/u)i yx'-/u)i yx^'/uji yx^"/u)i yx"/u) 

(5.74) 



which are obtained from (5.68), giving 



_, , _, ,X(a;28a)X(a;30a) 
^6(a)=^o(a' ^ ^ ^ ^ 




X{x^^a)X{x^^a) 



1 _ 2;60m.+26awi _ 60m+28awi _ 60m+32a\ 



X 



n — a;60™+2^)(l — a;60m+4aV]^ _ ^eOm+6 g.\ 
(I — 2;60m+34aW-| _ 60m+36aN/i _ 60m+38aN 






60m+42aNQ _ 60m+44a 



(1 _ x^0m+56a\fl _ rf.60m+58a\ I ' 



y6(l/a) = ^o(l/a): 



F(xiVa)F(xiVa) 




1 _ a;6'^™~3^ .^ V 1 — X^^'^~^^-Ml — a;60m-34b 
g/ '- a'^ a 

60m-14AVX _ 3.60m-12iV]^ _ 2;60m-8A' 
a' *> a'^ a' 

60m-28 6\/'1 ^60m-26b\ 



a' ^ a' ^ 



X 



*> a'^ a' ^ a' 

60m+16 6vi ™60m+18 6\"i 



^l_^60m+16i)(i 



X 



A 



6- 



b^ A(x2)A(xi2)A(a;20)A(a;30) 

X(x24u;)X(a;26u;)F(l/x26«;)F(l/x24^) 
X(a;6w)X(a;8w) F(l/a;8w)F(l/x6w) 



X 



30 



(5.75) 



Substituted into the eigenvalue expression 

w* A{x^w)A{x'^^w)A{x'^°/w)A{x^°/w) 



this gives (with 6 = — 1 and with elhptic nome x^°), 



(5.77) 



Ao 



E{-x^/w)E{-x'^lw)E{-x^'^/w)E{-x^^/w) 

E{-x^w),E{-x^w), E{-x^^w)E{-x^^w) 

E{-x^^w)E{-x^^w)E{-x^^w)E{-x^^w) 

E{-x^yw)E{-x^yw)E{-x^^/w)E{-x^^/wy 

At the isotropic point this reduces to 



■w 



X- 



Ae 
Ao 






■ 20 'J 



(5.75 



(5.79) 



5. 7 Mass m-j 



We begin with Wj = a^ for j = 1, . . . , A^ — 8 and wn^-j = bix , wn~6 = b2X , 



WN-n = 63a;" , WN-i = &4a; , wn-3 = 65a;' , wn-2 = bex^^ , wn-1 = brx~^, 
u!n = bgx^. We show in Appendix B.7 that the bi are all equal. The Bethe 
equation for the other roots is 



-co 






N 



g g/ga^ E{x%/a)E\x%/a)E\x^%/a) 
^ ^~^ ' ¥ E{x^a/b)E^{x^a/b)E^{x^^a/b) 



E^{x^%/a)E{x^*b/a) ^-^ E{x^'^a/aj)E{x^'^a,/a) 
^ E^{x^^a/b)E{x^^a/b) /J E{x^^aj/a)E{x^Oa/aj)' 

In the a; -^ limit this gives 

a^-8 + l(^^_g58)5/8/^8 ^ Q^ 



(5.80) 



(5.81) 



which is an (A^ — 8)th order equation for A^ — 8 zeros, so that there is again 
no hole. Equating this with YlfSi^{a — aj) we obtain 



-{A^-sb'f = A^.sb', 



(5.82) 



and from Bethe equations (B.40), 



UJ 



iA^-sbr^' 






l5N 



(5.83) 
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Defining 

7V-8 oo 

^rH=n n(l-^''V%), (5.84) 

j=l k=0 

Af-8 oo 

Gril/w)=l[l[il-x^^''aj/w), (5.85) 

j=i fc=i 

and rearranging the Bethe equation (5.80) as before, gives tlie recurrences 

^^""^ ~ A(x22a) X(x6a)X2(a;8a)X2(xi0a)X2(a;i2a)X(a:i4a) 

■^7(a:^°a) 
^J^7(a:20a)' 

_ A{x^^/a) Y{l/x'^a)Y^{l/x*a)Y^{l/x^a)Y\l/a)Y{x^/a) 
^^^/"^ ~ A(a;30/a)r(a;Va)F2(a;8/a)y2(a,io/a)y2(a;i2/a)y(a,i4/a) 

^7(a:^7a)- ^^ ^ 

Tlie solutions are 

X{x'^^a)X{x^^a)X{x^'^a)X{x^^a)X{x^^a) 



Tj{a) =J-'o{a 




X{x^a)X{x^Oa)X{x^^a)X{x^^a)X{x^^a) 

-)(1 — rF60m+28a\/'i _ „eOn 



X _ 2;60m+24a\Q _ 60m+28a\Q _ ^60m+30a^ 



_ x^Om+eaV^ _ a;60m+8aVX _ j.60m+10a\ 
X _ 2;60m+32aNQ _ ^60m+34a\Q _ ^60m+36a\Q _ ^60m+38a\ 



(5.87) 



-Vl — x60m+42aN|'i _ 60m+46a^ 



^X — 3;60m+60SiWx _ 2,60m+62a\Q _ ^60m+64a^ 



y7[ I ) yo{ / fY{x^ya)Y{x^ya)Y{x^ya)Y{x^ya)Y{x^ya) 

oo f Q _ a;60m-34bNQ _ ^60m-32fe\Q _ 60m-30b\ 

-'■-'- I fX — x^^"^^^^-)(l — a;^o™^i2Avx — x^Om-ioi) 

m=l liV g/V g/V g/ 

^ a'^ a'^ a'^ a' 



(1 - a;60'n-8.^)(l - a;60'n-6b)(X - x60™-4A)(X - a;60m-2b) 
^ a'^ g'^ g'^ a' 

'Y _ 2;60m+20 6Wx _ 2;60m+22 6\Q _ ^60m+24b 



CX — X^'^™-)(1 — X^'^'"+2AVX — a;60m+6A^ 



(5. 



32 



Substitution into 



Aj: 



-w°b- 



r^ .A{x^w)A{x'^'^w)A{x'^^/w)A{x^^/w) 



A{x^)A{x^^)A{x^o)A{x^o) 
X{x^'^w)X{x^''w)X{x^^w)X{x^^w)X{x^'^w) 
X{x'^w)X{x^w)X{x^Ow)X{x^^w)X{x^'^w) 
Y{l/x^^w)Y{l/x^^w)Y{l/x''^w)Y{l/x'"^w)Y{l/x'^w) 
Y{l/x^*w)Y{l/x^'^w)Y{l/x^Ow)Y{l/x^w)Y{l/x^w) 
xj^r{x^%)g7{l/x^%), (5.89) 



X 



X 



yields the result (with b = —1 and elliptic nome x 



60\ 



M_ ^ E{-x^/w)E{-x^/w)E{-x^^/w)E{-x^'^/w)E{-x^^/w) 
A^~^ E{-x*w)E{-x^w)E{-x^^w)E{-x^^w)E{-x^'^w) 
E{-x^^w)E{-x^^w)E{-x^^w)E{-x^^w)E{-x^^w) 
E{-x^^/w)E{-x^^/w)E{-x^^/w)E{~x^'^/w)E{-x^^/wy 

This reduces to 

2 



Ao 






60 



60 



(5.90) 



(5.91) 



at the isotropic point, with the elliptic functions of nome p 



8/15 



5.8 Mass rris 



We begin with the distribution Wj = aj for j = 1, 



A^ — 10 and wn-9 
b^x^^^ 



, Wn-4 



biX^^^, Wn-8 = b2X^^, Wn-7 = b3X^'^^, Wn-6 = b^x'^^, Wn-5 

bgx^^, Wn~3 = bjx^^, Wn~2 = bsx^, Wn-i = bgx^'^ , wn = biox'^ . We show in 
Appendix B.8 that the bi are equal. The Bethe equation for the other roots is 



-u 



Ejx^ya) 
'' E(x^^a) 



N 



rlo^5/8«^° Eix^b/a)E\x'^b/a)E^{x%/a) 
^ ^"'° ^ 610 E{x^a/b)E^{x''a/b)E'^{x^a/b) 



X 



E^{x^^b/a)E\x^%/a)E{x^^b/a) ^^V" E{x^'^a/aj)E{x^'^a,/a) 
E^{x^^a/b)E^{x^^a/b)E{x^^a/b) |J E{x^Oaj/a)E{x^^a/aj)' 

In the a; -^ limit this gives 



-,N-10 



+ HAN^iob'r^vb'' = 0, 



UJ 



(5.92) 



(5.93) 
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for the A^ — 10 zeros. We equate this with nj=i (a — aj) to obtain 

-{AM-iob^'f = A^^.ob'', (5.94) 

00 



and from Bethe equations (B.43), 






l6N 



fj6N 



1. (5.95) 



The functions 



N-IO oo 

FsM= n n(i-^''v«.), (5.96) 

j=l k=0 

N~10 oo 

Gsil/w)= n n(l-^'%/^)' (5-97) 

j=l fc=l 

must as a consequences of (5.92) obey the recurrences 



T^ia) 



Qs{l/a) 



A(a;22a) X{x^a)X'^{x'^a)X'^{x^a)X'^{x^^a)X'^{x^^a) 

X(x^^a) J^8(3:^°a) 
^X(a;i5a)J^8(a;20a)' 
A(a;^Va)y(l/a;^a)y'(l/x5a)y^(l/A) 
A(a;30/a) F(a;Va)F2(a;Va)F2(a;9/a) 
r2(i/a;a)F2(a;/a)y(a;3/^) g^^x^^a) 



Y^{x^^/a)Y\x^^/a)Y{x^^/a) Gs{x^Vo) ' 
Solving for them, we obtain 



._^, .X{x'^'^a)X{x'^^a)X{x^^a)X{x^^a)X{x^^a)X{x^'^a) 




X{x^a)X{x^a)X{x^^a)X{x^^a)X{x^^a)X{x^'^a) 

_ ^60m+25a\Q _ ^60m+27a\Q _ ^60m+29a\ 



X _ 2,60m+31a\Q _ ^60m+33aNQ _ ^60m+35a\ 






'\ _ 2;60m+35a\Q _ 60m+37a\Q _ ^60m+39a 



fX _ 2;60m+41a\Q _ ^60m+43aWx _ 2;60m+45a\ 
n — 2;60m+6laVX _ ^T^O^+^S^VX — a;60m+65a^ 



34 



(5.98) 



(5.99) 



(5.100) 



G^{l/a) = g^{l/a) 



Y{l/x^a)Y{l/x^a)Y{l/xa)Y{x/a)Y{x^/a)Y{x^/a) 
Y{x^^/a)Y{x^ya)Y{x^ya)Y{x^^/a)Y{x^^/a)Y{x^^/a) 




X 



60m-35 b 



-)(1 
a' ^ 



X 



60m-33 6 



-)(1 
a' ^ 



X 



60m-31 b-^ 



60m-15 b 
60m-29 b 



X'"""^^^-)il — X 
a' ^ 



-)(1 
a' ^ 



60m-13A)n — -r-eOm-ll b^ 

x—--'^)il-x- 

a' ^ 



(1 - a:60'»-9^)(l - a;60'n-76)(l _ a;60m-5 A) 
^ a' ^ a' ^ a' 



'I — X 



60m+15 b 



a' ^ 



X 



60m+17 b 



a' ^ 



X 



60m+19 b^ 



^ a' ^ a' ^ a' 



'I — X 



60m+21 b 



-)(1 
a' ^ 



X 



60m+23 b 



^Ml-x 

a' ^ 



60m+25 b ^ 



The eigenvalue expression is 



(5.101) 



A. 



w^ A{x'^w)A{x^'^w)A{x^^/w)A{x^^/w) 
Ifi A(a;2)A(a;i2)A(a;20)A(x30) 

X{x^^w)X{x^^w)X{x^^w)X{x'^^w)X{x^^w)X{x'^'^w) 
X{x^w)X{x''w)X{x^w)X{x^^w)X{x^3w)X{x^^w) 

Y{l/x^'^w)Y{l/x^^w)Y{l/x'^^w)Y{l/x'^^w)Y{l/x^^w) 



X 



X 



X 



Y{l/x^^w)Y{l/x^^w)Y{l/x^^w)Y{l/x^w)Y{l/x'^w) 
Y{l/x^'^w) 



Y{l/x^w) 



Tsix'''w)gsil/x'''w). 



(5.102) 



Thus the leading excitation in the w^ band, with b = —1 and elliptic nome 



x^°, is 



A 



8 _ QE{-xyw)E{-xyw)E{-xyw)E{-x^'^/w) 



X 



X 



E{-x^w)E{-x'^w)E{-x^w)E{-x^^w) 
E{-x^^ /w)E{-x^yw)E{-x^^w)E{-x^'^w) 
E{-x^^w)E{-x^^w)E{-x^yw)E{-x^yw) 
E{-x^^w)E{-x^^w)E{-x^^w)E{-x^^w) 
E{-x^^/w)E{-x*'^/w)E{-x^^/w)E{-x*yw) 



At the isotropic point this reduces to 



As 
Ao 



M-,)M^)Mfn)MfK)Mfu)Mo) 



M-,) M'-w) M¥) M'-w) Mtk) Ui 



H 2 



2' i 



(5.103) 



(5.104) 



where the elliptic functions are of nome p 
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6 General formula and correlation lengths 



Having obtained tlie relevant eigenvalues in the thermodynamic limit we are 
now in a position to calculate the correlation lengths and related mass gaps. 



6.1 Correlation lengths 



Recalling that our transfer matrix acts in the vertical direction, we consider 
the pair correlation function between two sites in the same column of the 
lattice separated by distance /. The correlation length defining the decay of the 
correlation function for large / can be obtained either by integrating over the 
relevant bands of eigenvalues and applying the method of steepest descent, or 
equivalently via the leading eigenvalue at the isotropic point (see, e.g., [19,20]). 
Here we simply follow the latter approach. 



Define the quantity 



rJu) = lim -^. 6.1 



The various correlation lengths follow as 

^7i = _logr,(M), (6.2) 

where we are to understand that we take the relevant leading eigenvalue at 
the isotropic point u = 3 A/2. 

6.2 Regime 1 



For L odd in regime 1, we derived the general result (4.18). For the leading 
eigenvalue in the band, our numerical checks confirm that the hole is located 
at a^ = —1, with the excitation parameter 6 = — 1, as expected. Thus using 
(4.20) gives the result 



r =2iog 



^4(ff,P"/'") 



(6.3) 



The correlation length diverges at criticality, with 



^^47!^"^' as p^O, (6.4) 
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where the correlation length exponent Vh is given by 

2(L + 1) 



bs 



3L 



(6.5) 



6.3 L = 3 regime 2 



The magnetic Ising value A = ||- for L = 3 in regime 2 has been of particular 
interest. Our results for the eight eigenvalue expressions obtained in Section 
5 can be summarised in the compact form 



r^.(u,)=u,"Wj] 



E{-x''w)E{-x^^+''/w)' 



(6.6) 



where the numbers a and n{a) are given in Table 2. The Eg numbers a have 
already appeared in [15] for the related Hamiltonian. The number n{a) de- 
notes the relevant band of eigenvalues. 

Table 2 



Parameters appearing in the eigenvalue expression (6.6). 



j 


n{a) 


a 


1 


2 


1, 11 


2 


2 


7, 13 


3 


3 


2, 10, 12 


4 


3 


6, 10, 14 


5 


4 


3, 9, 11, 13 


6 


4 


6, 8, 10, 12 


7 


5 


4, 8, 10, 12, 14 


8 


6 


5, 7, 9, 11, 13, 15 



In the original notation (6.6) reads 



rj{u) = n 



Q (an 8u „8/15\„Q ( air , 8u „8/15N 



"^3160 



||,p8/15)^^(|| 



8m 
15 



(6.7) 



where 'dz{u) = 'dilu + ^). The zero- momentum excitation energies of the 
related Hamiltonian, obtained by inverting the Fourier transforms in the ther- 
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modynamic Bethe Ansatz approach [11,15], are recovered on taking the loga- 
rithmic derivative of this expression evaluated at -u = and applying a Landen 
transformation. 



The first correlation length is given by 

-^,(l|,p8/15)^^(l|zr^p8/15) 



er' = 2iog 



^,(^^,psn5)^,(^^,ps/i5) 



(6. 



All eight fundamental correlation lengths can be written 



^. ^ ^-1 = 2 viog ^Q° ^ - feg) 



In particular. 



m^- ~8//^^^sin§ as p -^ . (6.10) 



This is the formula obtained by McCoy and Orrick [15], from which the E^ 
masses in (1.2) are recovered by virtue of trig identities. 

We see that there is surprisingly little variation in the masses (6.9) as a func- 
tion of the magnetic field-like variable p. In the limit p ^ 1 the mass ratios 
are given exactly by l|, 2, 2^, 3, 3|, 4, 5, which are to be compared with the 
Eg, values in (1.2). 



6.4 Universal magnetic Ising amplitude 



Making use of the Poisson summation formula in the free energy (3.19) at the 
magnetic Ising value A = || we find 

sm T 



/^4V3 ^-p^*^'^^ as p^O. (6.11) 



cos 30 



On the other hand, from (6.8) we have 

1 



ii- , fi . . P"'^'' as p^O. (6.12) 

8v3 sm| 
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Combining these results gives the universal magnetic Ising amplitude 

/C? = j= = 0.061728 589... as p^O. (6.13) 

16V<jSm g cos gQ 



This result has been predicted earlier by other means. Namely by thermody- 
namic Bethe Ansatz calculations based on the E^ scattering theory[26,5,7,27] 
(see also Ref. [28] in the context of the form-factor bootstrap approach). Here 
it is obtained explicitly from the lattice model. 



7 Conclusion 



We have applied the exact perturbation approach to the Bethe Ansatz solu- 
tion of the dilute Al lattice model to derive the free energy per site and the 
excitation energies in regimes 1 and 2 for L odd. In the dilute Al model the 
elliptic nome p is magnetic field-like for L odd and is temperature-like for L 
even. The particular point in regime 2 for L = 3 has attracted considerable 
recent attention, being in the same universality class as the Ising model in 
a magnetic field. We have specifically considered the case L odd for which 
the method perturbs from the ordered high field limit. Our result for the free 
energy (3.19) is in agreement with that obtained via the inversion relation 
method for general L [8,9]. 

In regime 1 the leading excitations in the w band are 1-strings, characterised 
by the string excitation parameter h and a single hole. Our final result for the 
related correlation length is given in (6.3). The excitations are considerably 
more complicated in regime 2, where we have concentrated on the case L = 3. 
Here our results for the leading eigenvalue associated with each mass are 
summarised in (6.6) and (6.7). The corresponding inverse correlation lengths 
or masses are given in (6.9), all of which give the magnetic Ising correlation 
length exponent z//i = ^. In particular, the E^ masses (1.2) predicted by 
Zamolodchikov appear with the approach to criticality. 

Of course these masses have been obtained earlier in the scaling limit via 
the thermodynamic Bethe Ansatz approach [11]. Although the calculations 
are somewhat complicated our approach nevertheless provides in principle 
a means of classifying and counting all excitations in the eigenvalue bands. 
We have not pursued this classification in detail, rather contenting ourselves 
with the location of the leading eigenvalue in each band, of relevance to the 
correlation lengths. As for the thermodynamic Bethe Ansatz, our calculations 
rely upon the key input of the string configurations. We have taken as our 
starting point the eight thermodynamically significant string excitation types 
revealed in the previous numerical studies [11,16,17]. The one exception is that 
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we have considered the 7-string configuration [17] for mass m^ (see Table 1) 
rather than the 5-string [11]. However, although we have not included here our 
calculations based on the 5-string, we find that the 5-string configuration is 
also consistent in the ordered limit. Moreover, we find that both the 5-string 
and the 7-string lead to the same eigenvalue expression and thus mass. The 
origin of this behaviour is seen in the 7-string calculation where the string 
parameters a and /3 (see eqn (B.30) do not appear in the final result (5.54). 
In both cases, it is the unpaired string at bx^^ which ultimately drives the 
calculation. 

The derivation of the correlation length for L 7^ 3 in regime 2 is complicated. 
In this regime the leading excitation in the w band has £ = |(L — 1) + 1 and, 
like the leading 2-string in the w"^ band for L = 3, it begins life for small N and 
p ~ as a 1-string. We have not pursued this further. Nevertheless we have 
numerically observed that the final result (4.19) for regime 1 also applies to 
the leading w band excitation in regime 2. We thus believe that the correlation 
length (6.8) and the corresponding exponents 

r 2(L + 1) , , 



also hold in regime 2 for L ^ 3. This result includes Uh = j^ for L = 3. 

The correlation length exponents are seen to satisfy the general scaling relation 
2z//i = 1 + 1/6, which follows from the general relation /^^ ~ constant, where 
/ ~ pi+i/^ is the singular part of the bulk free energy and the exponents 
6 are those following from the singular behaviour of (3.19) [8,9]. The same 
correlation length exponents should hold for L even, for which the integrable 
perturbation is thermal-like. The scaling relation is now 2z/i = 2 — a, where 
Ut is as given in (6.5), (7.1) and a in [8,9]. In particular, (6.5) gives the Ising 
value Vt = 1 for L = 2 in regime 1, as expected. 

Another approach to the results given here is to use the inversion relation 
method, as was originally used to obtain the bulk [8,9] and excess surface [10] 
free energies. This approach can also be applied to the excitations. This has 
been done, for example, for the eight-vertex [30] and the Andrews-Baxter- 
Forrester models [31]. Although in principle completely avoiding the string 
hypothesis the inversion relation method still requires strong assumptions on 
the analyticity properties of the eigenvalues. Our explicit results for rj{u) 
give these analyticity properties a posteriori. It remains to be seen if these 
properties can be established a priori, thus allowing an alternative derivation 



"^ Yet another would involve integral equations for the root densities, as done for 
the eight- vertex model [29]. 
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of the mass gaps. We note that the inversion relation 

rj{u)rj{u + 3\) = 1 (7.2) 

is indeed satisfied by our results. Alternatively, rj{w) rj{x^^w) = 1 for L = 3 
in regime 2, which is seen to hold trivially in view of (6.6). There is a further 
relation^ 

r,{u)r,{u + 2X)=r,{u + X), (7.3) 

which is also easily seen to be satisfied by our results. 
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Appendices 

A Auxiliary functions 

The simplest identities for the auxiliary functions (4.6)-(4.8) are 



1-xw/b), (A.l) 



F(x"/m;) 
R{x'^w) 

5(x"-27w) 



(l-a;"6/w;), (A.2) 

(l-x^w/oTv), (A.3) 

(l-x"a7v/w). (A.4) 



SiyX'^/w) 
When solving the recurrence relations (in either regime), their generalisations 



oo T^/ 12ms+n \ oo 

n w;L..„J.L = n (1 - ^-'^'^-yn), (a^s) 

m=0 ^ \-^ ^ ) m=0 

oo -y ( ^12ms+n—2r l„.,\ oo 

n |l.a„..../:r' = n (1 - x""-»v-), (A.6) 

m=0 -' l-^ / "^/ m=0 



are used. In regime 1 we also require 



oo JD(r^l2ms+n .\ oo 

n7jP;;3T;s&=n(i--""''"'Va.). (a.t) 

m=0 -"-l-^ "^y' m=0 

oo c/ 12ms+n-2r /„,,\ oo 

m=0 "^l^ /^J «=0 

In regime 2 we treat only L = 3, so that s = 5 and r = 16. For mass m^ we 
need the analogous identities, 



oo p (^mm+n„ ,\ oo 

n fil'L.n.3"! = n (1 - ^»"'^"»/«A'-C,). (A^9) 

oo c /™60m+n-32 /„,,\ oo 

n 5l«o..../ii n(i - .""■-"<.«-«/-). (A,io) 

Of course, the regime 2 auxiliary functions X^iw) and X^(u') obey the same 
identities as X(w), with h therein replaced by a and /9, respectively. Simi- 
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larly Yo,{l/'w) and Yi3{l/w) obey the same identities as Y{l/w), with again h 
replaced by a and {3 as appropriate. 

When substituting JF and Q into the k = N Bethe equation in regime 1 we 
use 



oo -i^/ (2r-4s)«+n-2r/ \ oo 

oo p/'„(2r-4s)«+n \ oo 

oo q( {2r-4s)l+n-2r /„,,\ oo 



To combine the Bethe equations from the strings in Appendix B, once the 
functions jFj and Qi have been found, we use 



oo A(^60m+n/l^\ oo 

n 2,L.J:L = n (1 - -"'-'V")". (a.i5) 

m=0 ^y-^ I") m=0 

oo A(^QGm+n-L\ oo 

n 2L....X - n (1 - -'^"'""r, (a.w) 

m=0 ^l-^ "l m=0 

oo 4 / 60m+60+n //,\ oo 

n %.,.,.JZ -x(-'r'' n (1 - x«""*>«*" lb)" , (A,i7) 

m=0 ^l-^ /"/' m=0 

oo 4 /„60m+60+nl,\ oo 

n i&ss^ =>'(-")-" nil--™"*"*"'')" (A.18) 

m=0 ^l-^ "^ m=0 

These last four relations are just special cases of more general identities (which 
we do not require) which would have 12s in place of 60 and 2r in place of 32. 



B Bethe equations for the masses 



In this Appendix we consider the Bethe equations for the eight string types 
corresponding to the E^i, masses. Their chief usefulness is to establish ana- 
lytically that the coefficient of each member of a string is the same. (In the 
case of m4 and me, these equations give equahty only between subsets of the 
coefficients, and this feature was explicitly included in the calculations in the 
body of the paper.) Secondly, once the auxiliary functions (which give the re- 
mainder of the Bethe roots, the a^) are substituted into the string equations. 
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it is possible to establish a condition on the string coefficient b in which the 
same patterns in the exponents of x arise as do in the eigenvalues (or masses) 
themselves. Such higher level Bethe equations for the string parameters have 
been discussed, for example, for the CSOS model [20]. In each case these equa- 
tions ensure that the corresponding eigenvalue expression reduces to an A^th 
root of unity at w = 1. This must be so, as the row transfer matrix reduces 
to a shift operator at this point. 

Because some of the strings are quite long, it will be convenient to introduce 
the notation ni^Li h = Bm- 



B.l Mass mi 



With the roots Wj = aj ior j = 1, . . . , N — 2 and wn-i = hix , wn = b2X 
the Bethe equations for k = N — 1 and k = N are 



11 



-U! 



X 



-10 



E{x^'/h] 
E{x/h) 



N 



iA,,^,B,)-'/%^x- 



N-18 



X 



E{x^%2/bi)E{x^%/b 



N-2 

n 



E{xaj/bi)E{x^-^aj/b 



E{bi/b2)E{x%2/bi) f^\ E{x^^a^/bi)E{x^3a,/b,] 



(B.l) 



—UJ 



X 



10 



E{xb2 



N 



E{x^^b2 
E{bi/b2)E{x%2/bi 



{A^^2B2)-''X^''^'' 



^^ E{x^^b2/a,)E{x^a,/b2) 

" E{x^%2lbi)E{x^%2lbi) /J E{xaj/b2)E{x^^aj/b2) ' 

Taken individually as a; ^ 0, these equations imply that 



;i - ^) = o{x'^) 



02 



fB.2) 



(B.3) 



Thus in the thermodynamic {N -^ oo) limit 6i = 62 = b. 

Forming the product of (B.l) and (B.2), we see that many factors cancel and 
those remaining can be written in terms of the auxiliary functions (5.5) and 
(5.6) for mi, namely 



E{xb)E{x^^/b) 
E{x/b)E{x^^b) 



N 



^2N 



j^i{x%)j^i{x^^b)gi{i/x%)gi{i/x^^b)' 



(B.4) 
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Using the functions (5.9) and (5.10) and identities from Appendix A, we obtain 



E{xb, x^^)E{x^\ x^^)E{x^^/b, x^^)E{x^^/b, x^^) 
E{x/b, x60)E(xiV6, x^^)E{x^^b, x^^)E{x^^h, x^o) 



-\N 



l2N 



(B.5) 



This defining relationship is the higher level Bethe equation to be satisfied by 
the string parameter b. Note that it is trivially satisfied by 6 = — 1. Compare 
also the pattern of exponents in this equation with the expression (5.12) for 
Ai. Taken together they ensure that the eigenvalue is an A^th root of unity at 
w = 1. 



B.2 Mass 1712 



With the roots tfj^ = a^ for j = 1, . . . , A^ — 4 and wn^^ = bix~^, wn_2 = b2X^ , 
wjy-i = b^x'^^, wjy = b^x^^ the Bethe equations for k = N — 3, N — 2, N — 1, N, 
respectively, are 



—u 



bix 



_^E{x^ybi) 
Eix%i) 



N 



{AM-,B,f'%lb, 



-^x-^^ 



X 



E{bi/b2)E{x^%2/bi)E{x^%^/bi)E{x'%i/b^)E{x%i/bi)E{x'%i/b 
E{x^%i/b2)E{x^%2/bi)E{b3/bi)E{x^%i/b3)E{bi/b^)E{x^%Jbi] 

^^ E{x%/aj)E{x^^aj/bi) 
^ 11 



t^ E{x^^aj/bi)E{x^%i/aj 



:b.6) 



-LO 



b2X 



,E{x'/h) 



-\ N 



Eix^^bo 



{Aiy^4B4 



,2^-2^18 



2'^3 



X 



X 



E{x'%,/b2)E{x^%2/b,)E{b;,/b2)E{x^%2/b3)E{b2/b^)E{x^%Jb2) 
E{bi/b2)E{x^%2/bi)E{x%2/b3)E{x^%2/b3)E{x^%2/b4)E{x^%Jb2, 
^/^'' E{x'%/a,)E{x^^a,/b2) ^ 



r- 



^^ E{x^aj/b2)E{x^^b2/a 



Jy 



B.7) 



-u 



X 



-10 



Eix^^bs) 



N 



iAj,_,B,f^A],'_,b, 



2N-6 



bfx-'''-' 



X 



Eixyb^) 

E{bs/bi)E{x^%Jb3)E{x%2/b3)E{x^%2/bs)E{x^%4/b5)E{x^%Jbs) 



E{x^%3/bi)E{x^%i/bs)E{bs/b2)E{x^%2/b3)E{x%4/b3)E{x^%4/b3 
^^ E{x^^aj/bs)E{x^aj/b3) ^ 



^\ E{x^^aj/b3)E{x%/a 



B. 



'3 J 
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—u 



X 



10 



E{x%) 
E{x^) 



N 



{Ar,-,B,f/'A-^'_,br-%^'x^ 



2N-67-2HN+6 



X 



E{bi/h)E{x^%i/bi)E{x^%2/h)E{x'%i/b2)E{x%i/b3)E{x^%i/h 
E{x%i/bi)E{x^%i/bi)E{b2/b4)E{x^%Jb2)E{x^%4/b3)E{x^%i/b3 

^^ E{x^aj/bi)E{x^aj/b4) 
^ f}^ E{x^^aj/b^)E{x%Jaj)' ' 

Taken individually as x ^ 0, these equations imply that 



B.9) 



(1 


-'m- 




= 0(x^^), 


(1- 


bz\ 
b2> 




(1-&) 




(1- 


bj 


(1 




bj 


= 0(a:^^), 


(1- 
(1- 


b4) 

bj 



Oix 



2N\ 



Oix 



2N\ 



(B.IO) 



We cannot read off the relationship between the 6j in this case, as we could 
for mi, but we see that these relations are satisfied if 



3N\ 



&i = 64 + 0(x^^'), 61 = 62 + 0(0;^^), bi = b3 + 0ix 



„3Af^ 



(B.ll) 



Thus we conclude bi = b2 = b^ = b^ = b and 84^ = b'^. 

The product of the four Bethe equations can be written in terms of the aux- 
iliary functions (5.19) for 1712-, with 



E{x^yb)E{x^yb) 



1 N 



E{x^^b)E{x^^b) 



2N 



T2ix'b)T2ix%)T2ix'%)T2ix"b) 



-1 



Q2ii/x'b)g2ii/x%)g2ii/x''b)g2ii/x''b) 



-1 



[B.12) 



Using the functions (5.21) and (5.22) and identities from Appendix A, we 
obtain 



Ejx'^b, x^'')E{x^%, x'^'>)E{x^yb, a;60)E(a;^V&, x^") 
E{xyb, x^o)E{x^yb, a:60)E(x376, x(^^)E{x^%, x^^) 



N 



u2N 



[B.13) 



B.3 Massms 



With the roots Wj = Uj for j = 1, . . . , A^ — 4 and 10^-3 = bix , wn-2 = b2X^, 
wn-1 = &3a;~^°, wn = &4a;^° the Bethe equations for /c = A^ — 3, A^ — 2, A^— 1, A^ 
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read 



—LO 



h\x 



-10 



E{h) 



N 



{AM-,B,f'%lh 



-2^-38 



X 



E{x%2/bi)E{x^%2/bi)E{x^%3/bi)E{x^%i/b3)E{x^%Jbi)E{x'^%Jbi 



E{bi/b2)E{x^%2/bi)E{b3/bi)E{x^%i/b3)E{x%4/bi)E{x^%i/bi 
^-^ E{bi/a,)E{x^'>a,/b,' 



t^ E{x^Oaj/bi)E{x^%i/aj] 



B.14) 



—u 



X 



X 



10 



E{b2 



Eix'^%, 



N 



{AN^^B,f'%lbl 



-2^38 



X 



E{b,/b2)E{x^%2/bi)E{x%2/b3)E{x'%2/b3)E{b2/b^)E{x'%Jb2) 



E{x%2/bi)E{x^%2/bi)E{x^%2/b3)E{x^%2/b3)E{x^%2/b4)E{x^%Jb2) 
't^^E{x^%2/ai)E{x^^a,/b2) 



i=i 



E{aj/b2)E{x^%2/aj 



[B.15) 



— CJ 



X 



a; 



AT 



iAr,^,B,)'/'A],\br-'x~''''^'' 



Eix^ybs) 

E{b3/bi)E{x^%i/b3)E{x^%2/b3)E{x^%2/b3)E{x'^%4/b3)E{x^%i 



E{x^%3/bi)E{x^%i/b3)E{x%2/b3)E{x^%2/b3)E{x^%4/b3)E{x^%Jb3) 
^^ E{x^''a,/b3)E{x^aj/b3) 



f^, E{x^^a,/b3)E{b3la,) 



(B.16) 



-u; 



X 



10 



E(a;3064) 



AT 



{AM-.,B,f"'A-^\b 



27V-8^18Af-32 



X 



E{x%4/bi)E{x^%4/bi)E{x^%2/b4)E{x^%i/b2)E{x^%4/b3)E{x'^%4/b3 
E{x^%Jbi)E{x^%Jbi)E{b2/b4)E{x^%Jb2)E{x^%Jb3)E{x^%Jb3) 
^-^ E{x''%4/baj)E{aj/b4) 



^ jii E{x%Jaj)E{x^%Jaj)' 



B.17) 



Taken individually as x -^ 0, these equations imply that 



&2^^ b,' 



bi 



)(1 



62. 






:-|) 
:-|) 






fB.18) 
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from which we can read off 

6i = 63 + 0(x^^), 62 = &4 + 0(x^^) 



61 = 62 + 0(x 



2N\ 



(B.19) 



Thus we conclude bi = b2 = b^ = b^ = b and B4 = b'^. 

Expressing the product of the four Bethe equations in terms of the rris auxihary 
functions (5.30) gives 



E{x^yb)E{x^yb)E{x^yb) 
E{x^%)E{x^%)E{x^%) 



1 N 



l3N 



J^^{x%):Fl{x^%):F,{x^'b)G,{l/x%)G'^{llx^%)G,{llx^'b) . (B.20) 



Using the functions (5.32) and (5.33) and identities from Appendix A we 
obtain, in terms of elhptic functions (of nome x^°, not x^^), 



E{x'b)E{x'%)E{x^^b)E{x'yb)E{x^yb)E{x^yb) 
E{xyb)E{x^yb)E{x^yb)E{x^%)E{x^%)E{x^^) 



N 



U3N 



(B.21) 



B.4 Mass 1714 



With the roots Wj = aj ior j = 1, . . . , N — 7 and wn^g = bix \ wtv-s = b2X, 



wn-4 = hx ^, wn~3 = b^x'^, wn-2 = hx 
Bethe equations ior k = N — 6, . . . ,N are 



-13 



w^^i = bgx^^, wn- = bjx^^ the 



-LO 



bix 



.iE(a;"/6i) 



N 



{AM-7Bjf'%\{b,b,b-j)-^x 



-2^-6 



X 



E{x%i) 

E{x%i/b2)E{x^%i/b2)E{x^%/bi)E{x^%i/b^)E{x%i/bi)E{x'^%4/bi) 
E{x^%2/bi)E{x^%i/b2)E{x%/bi)E{x^%i/b3)E{x^%4/bi)E{x^%4/bi] 
E{x%/bi)E{x^%5/bi)E{x%6/bi)E{x%6/bi)E{x'^b7/bi)E{x%/bi) 



E{b^/b^)E{x%i/b^)E{x%i/bQ)E{x^bi/b(;)E{x^bi/bj)E{x%i/bj) 
y E{x%ilaj)E{x^^ajlb^) 
^ y E(xiia,/6i)E(xi%/a,)' ' 



B.22) 



—OJ 



b2X 



E{x^/b2 
Eix^^b. 



N 



{AN--jB,f'%l{b,b,h 



-2x6 



E{x^%2lhi)E{x^%i/b2)E{x^%2lh)E{x^%^lb2)E{x%/b4)E{x'^%Jb2) 
E{x%i/b2)E{x^%Jb2)E{x%3/b2)E{x^%2/b3)E{x^%Jb2)E{x^%Jb2) 
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E{x%/b5)E{x%5/b2)E{x%/b2)E{b2/be)E{x%/b2)E{x%/b2) 

^E{x%2/b5)E{x%2/b5)E{x%2/b6)E{x^%2/bQ)E{x^2/b7)E{x%2/b7 

^-^ E{x^%/aj)E{x^^aj/b2^ 
^ 11 



^\ E{x^a,/b2)E{x^^b2/ 



B.23) 



'jy 



-00 



63X" 



■Eix^-^/b^) 



N 



{Aj,^rBjf%tib,b,)-'x- 



34 



X 



X 



E{X%3) 

E{x%/bi)E{x^%i/b3)E{x%3/b2)E{x^%2/b3)E{x%/b3)E{x%4/b3) 
E{x^%3/bi)E{x^%i/b3)E{x^%2/h3)E{x^%3/b2)E{x%3/hi)E{x%3lh) 
E{x'%/b3)E{x'%/bs)E{x^%,/b3)E{x%,/b3)E{x'%/b3)E{x'%/b3] 



E{x%,/b,)E{x%,/b3)E{x%3/k)E{b3/b,)E{xbr/b3)E{x%r/b3 
^^ E{x%/a,)E{x^''aj/b3' 



^^ E{x^^aj/b3)E{x^%3/aj 



B.24) 



-UJ 



64a; 



-jEix^b^) 



-\ N 



iAN^rBrf%lib5bs)-'x 



34 



X 



X 



E{x^'^bi) 

E{x^%Jbi)E{x^%Jbi)E{x^%Jb2)E{x^%Jb2)E{x%3/b4)E{x%3/b4) 
E{x%i/bi)E{x^%i/b^)E{x%2/W)E{x^%Jb2)E{x%Jb3)E{x%Jb3) 

E{x%5/h)E{b5/b4)E{x%/bQ)E{x%4/b6)E{xb4/br)E{x%/br) 
E{x^%Jb5)E{x%Jb5)E{x^%JbQ)E{x^%JbQ)E{x^^bJbr)E{x^%Jbj 

""-' E{x^''bJaj)E{x^^aj/bi) 



n 



\ E{x^aj/b4)E{x'^'^b4/aj 



(B.25) 



-LU 



X 



-10 



Eix^^b^) 



N 



{A^^rBrf'A-^'^.b, 



2N-8 



{b,b2b,r'x-''-''' 



Eix^b,] 

E{b,/b,)E{x%,/b,)E{x%/b,)E{x%2/h)E{x%/b,)E{x%,/b3) 
E{x%5/bi)E{x^%5/bi)E{x%2/b5)E{x%5/b2)E{x^%/b3)E{x^%r^/bs) 
E{x^%4/b5)E{x%i/b5)E{x^%Q/b5)E{x^%/b5)E{x^%r/b5)E{x^'^b7/b 

E{x^,/b4)E{b5/b4)E{x%,/b,)E{x%,/b,)E{x^r/b5)E{x%j/b5) 
^-^ E{x^aj/b5)E{xaj/b5'' 



n 



\ E{x'^^aj/b^)E{x'^b^/aj 



B.26) 



—00 



X 



10 



E{x%q) 
E(x^%, 



N 



iAr,-jBrf'A-^'_,b, 



2N-8 
6 



{b,b2b,)-'x''+''' 



X 



E{x%,/b,)E{x%,/b,)E{x%2/be)E{x'%2/be)E{x%/b,)E{b3/be) 
Eix%,/b,)Eix%,/b,)Eix%,/b2)Eib2/b,)Eix^%,/b,)Eix%/b,) 
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X 



E{x^%JbG)E{x^%/be)E{x\/b5)E{x%G/b5)E{x'^bQ/bj)E{x^%/be) 



E{x%4/b6)E{x%4/b6)E{x^%6/b5)E{x^%e/b5)E{x^%7/be)E{x^'^be/b7 

^^ E{x^%/aj)E{x''aj/bG' 
^ 11 



j=i 



E{x%Q / aj)E{xbQ / a 



B.27) 



1) 



-UJ 



X 



E{X%7 



-\N 



{A^-,B,f''Af_,bf-'\b,b, 



E{x^%7) 
E{x%i/bj)E{x%i/h)E{x'^b2/h)E{x%2/b7)E{xbj/b^)E{x%j/b^) 



X 



E{x'^bT/bi)E{x%/bi)E{x%/b2)E{x%7/b2)E{x^^b7/bs)E{x^^b7/b3) 
E{x^^b4/b7)E{x^%i/b7)E{x^r/b5)E{x%r/h)E{x^%7/be)E{x^^bQ/b7) 



E{xb4/b7)E{x%Jb7)E{x^%7/b^)E{x^'^b7lb^)E{x'^b(^/b7)E{x^%7/bQ 

^^^^^^ E{x^%7/ai)E{x'a,/h7) 
f}^ E{x%7/aj)E{x%/aj) ' 

Taken individually as a; ^ these equations imply 



B.28) 



:i-^)=0(x^), 



(1 









06 









(B.29) 



which fall into two unconnected groups. Notice also that the seventh equation 
(B.28) does not provide any link between 67 and the other 6j. We are only able 
to conclude that 



bi 



a, 



^3 = ^6 = P, and 67 



(B.30) 



with B7 = a-^P% 



B.5 Mass m^ 



We consider the roots Wj = aj for j = 1, ... ,A^ — 6 and wat-s = bix~^^, 
W7V-4 = ^23^^^! "U^Af-s = ^32;"^^, W]\f^2 = b^x^^, Wj\f^i = b^x^^, Wj^f = b^x^ . In 
the X ^ limit, the Bethe equations for the latter six roots are 



X 



-1 



-^°^ - —{Aj,^,B,f''A-^\br-''x-'''-^'E{b„ 
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X- = —{A^^eB,f'A],'_,br~V^''E{h/b,), 



aoAf 



X 



-ION 



-1 

UJ 



{AM-,B,f"'A-^\h^-%fx-''-''^E{h/h)E{h,/h)] 



-1 



X— = —iAN^eB,f'A],'_,b^-%fx''-^''Eih/h)Eih/h), 



hx 
bex^ 



ION _ '_ 

UJ 
9 ^ -1 



— iAN^eB,f%t{b,be)-'x-'' [E{b,/b,)E{b2/b,)r' , 

UJ 

"" = —{Aj,^,B,fHl{b,b,)-'x'^E{b,/b,)E{b,/b,). 

UJ 



We conclude from the first, third and sixth of these that 

6i = 66 + 0(a:6^), b^ = b^ + O^x^""), 63 = 64 + 0(0;^^) 



(B.31) 



(B.32) 



From the second and fifth we have 



&2 = &5 + 0(a;^'"), b^ = b^ + 0{x' 



3N 



(B.33) 



The fourth confirms 63 = 64 + 0{x^'^), so that we may conclude that the bi 
are equal with Bq = b^. 

In terms of the rrir, auxiliary functions (5.63) and (5.64) for the product of the 
six Bethe equations give 



E{x'yb)E{x^^/b)E{x^yb)E{x^yb) 

E{x^%)E{x^^b)E{x^%)E{x^%) 
J^,,{x^b)J^l{x%)J^l{x^H)J^^{x^^b) 



N 



^AN 
X 



-1 



g,{i/x'b)gl{i/x%)gl{i/x''b)g,{i/x''h) 

The parameter b is thus found to obey the equation 



(B.34) 



N 



E{x%)E{x%)E{x^^b)E{x^%) 
E{xyb)E{xyb)E{x^^/b)E{x^yb) 

E{x^yb)E{x^yb)E{x^yb)E{x^yb) 
E{x^%)E{x^%)E{x^^b)E{x^%) 

where the elliptic functions are of nome x^^. 



N 



uiN 



(B.35) 
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B. 6 Mass rriQ 



We consider the roots Wj = a^ for j = 1, . . . , A^ — 7 and w^^q = bix , 
wn-5 = b2X^^, wn-4 = hx''^^, wn-3 = b^x^^ , wn-2 = hx"^, Wn^I = bQx'^, 
wn = bjx^^. In the x ^ hmit, the Bethe equations for the latter seven roots 
are 

x-''"" = —{A^^,Brr/'A-^%br-'%fx~'^~''E{b,/b,)-\ 



UJ 

X- = —iAr,^rBr)'/'A],'_,br-'%fx'''^''E{b,/b,), 



JON 



b,x-'' 
b^x^ 



"" = —{A^_,Br)'/%{bAr'x-''E{bs/b,r\ 

UJ 



N 



-1 



UJ 
N _ -1 

UJ 
N _-l 

UJ 



{AN--jBjf'%\{b,b,r^x'''E{b,/b,), 



{Ar,.yB^f%l{b2b,b,y 
{AN^jBjfl%l{b^b,b^y 



^x~'' 



Eibe/b, 



E{b,/b5)E{b,/b2y 
_2 ^oE{bi/be)E{be/bj) 

Jb 



E{b,lb, 



WN ~^ ( A n \5/S^WN^ih/h) 

UJ E{bQ bj) 



(B.36) 



We notice immediately that the third and fourth equations decouple from the 
others to give 63 = 64 + 0(x^°^), so we set 63 = 64 = «• From the remaining 
five equations we may conclude that (to 0(x^^) or 0(x^^)) the remaining bi 
are equal, with B-j = b^a^. 

The product of the seven Bethe equations may be written in terms of the 
auxiliary functions (5.72) to give 



X 



X 



E(x^7&)E(x=^V6)E(x^7a) 
E{x'^^b)E{x'^%)E{x'^^a) ^ 
'F6{x^%)Fe{x^%)Ee{x^%)Fe{x^%)Ee{x^%) 

FQ{x%)FQ{x%)Fe{x%)FQ{x%)Fe{x^%) 
\{l/x^%)Ge{l/x^%)Ge{l/x^%)Ge{l/x^%)Ge{l/x'%) 



^47V 



G^{\l x^%)Gq{\I x%)G^{\l x%)G^{\l x%)G^{\l x'^b) 

FG{x^a)Fe{x^^a)FG{x^^a) Ge{l/x^a)Ge{l/x^^a)Ge{l/x^^a) 
F6(xi0a)F6(xi2a)F6(a;30a)G'6(l/a;i0a)G'6(l/a;i2a)G'6(l/a;30a) 



Now, while some of these factors may be grouped to give the functions 
Qq in (5.73), for the others we must use 



(B.37) 
jFg and 
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Fe{w) = llJ^,{x''^w] 



XJx'^'+^ 



w) 



1=0 



X„fa;i2'+io 



w] 



oo Y / 121+2 I \ 



(B.38) 



/=o 



which is similar to what we did in regime 1. By applying identities from 
Appendix A, we find that all factors involving a cancel. The final result, in 
terms of elliptic functions of nome x^^ ^ is 



N 



E{x%)E{x%)E{x^%)E{x^%) 
E{xyb)E{xyb)E{x^yb)E{x^yb) 

'E(a;36/&)E(a;38/&)E(a:^V6)E(x^76) 
E{x^%)E{x^%)E{x^^)E{x^%) 



N 



L4:N 



[B.39) 



B. 7 Mass mj 



For the roots Wj = aj ior j = 1, . . . , N — 8 and Wn^j = bix , Wn-6 
wn-5 = bsx'^^, wn-4. = M^^ wn~3 = b^x'^^ , wn~2 = bex^^, wn^i -- 
u!n = bsx^ the last eight Bethe equations in the x — *> limit give 



-- b2X^^, 

bjx~^, 



X 



^lON 



AON 



-1 



— iA,,.sBsf^'Aj,\bl'''''x-'''-''E{bs/b^r\ 



-1 



X- = —{Ar,^sBsf''A-^\bf-'^x'^^^'E{b,/b,), 



X 



-ION 



UJ 

-1 

UJ 



{Ar,-,B,f"A~^\br-'%, 



2 -2N-Q2 



[E{b^/b,)E{b,/b^)]-\ 



X 

b,x-'' 
bsx^ 



ION 



iA,,^sBsr^'A-^\br-'%'x'''^''E{b,/b,)E{br/b,), 



UJ 

^ = —{Ar,-sBs)'^%{b,bs)-'x~'' [E{b,/b,)E{b,/b,)]-' , 

UJ 
N -1 



UJ 
N _ -1 

UJ 
N -1 



{AN-sBs)'^X{b5bj)-'x''^E{be/bs)E{b,/b,), 
{A^^sBs)'/mb,b,bsr'x''' [E{br/b,)E{b2/br)]-' , 
— iAr,^sBsT^Xihb5b7r'x''Eibs/b,)Eib,/bs). (B.40) 



UJ 



r.6N\ 



In this case we can conclude that the bi are equal, to 0{x ) or 0{x' 

The product of the eight Bethe equations can be written in terms of the rrij 
auxiliary functions (5.87) and (5.88) as 
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E{x^yb)E{x^yb)E{x^yb)E{x^yb)E{x^^/b) 
E{x^%)E{x^%)E{x^^b)E{x^%)E{x^%) 



AT 
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J^7{x%)J^^{x%)J^^{x'%)J^^ix'%)J^7{x'%) 



X gj{i/x%)g^{i/x%)g^{i/x'%)g^{i/x^^b)gr{i/x^%) 

Or finally, again with elliptic nome x®°, 



(B.41) 



N 



E{x%)E{x%)E{x^%)E{x^%)E{x^% 
E{xyb)E{xyb)E{x^yb)E{x^yb)E{x'^yb) 

E{x^yb)E{x^yb)E{x^'>/b)E{x^yb)E{x^yb) 
E{x^'^b)E{x^%)E{x^%)E{x^%)E{x^%) 



iN 



l5N 



(B.42) 



B.8 Mass ms 



For the roots Wj = aj for j = 1, . . . , A^ — 10 and Wn-9 = bix , w^^s = b2X , 

Wn~7 = bsX'^^, Wn-6 = b4^X^^, Wn~5 = b5X'^'^ , Wn-4. = b^X^'^ , Wn-3 = M~^ 

'Wn~2 = bsx^, wn~i = bgx"'^ , wn = biox"^ the last ten Bethe equations in the 
X ^ limit ffive 



X 



-2N 



-{Am- 



J2N 



X 



X 



-6N 



X 



6N 



X 



-9N 



X 



9N 



b-jX~ 



bgx' 



bwx'^ 



N 



N 



N 



N 



— {Am- 

— {An- 

UJ 

— {An- 

— {An- 

UJ 

— {An- 

UJ 

— {An- 

UJ 

— {An- 

UJ 

— {An- 

UJ 

— {An- 



lo-Bio 
lo-Bio 

lO-SlO 
lO-Sio 
lO-SlO 

lo-Bio 
lo-Bio 
lo-Bio 

lO-Sio 

lo-Bio 



'''A-^\,br~^'x-^'E{brolb,)-\ 
'/'A],\,br-VE{b,/b,), 
'"A-^\,br-^%-,^x-^>' [E{b,/b,,)E{b,/b,)]-' , 
"'A-^\,br-'%fx^^E{b,/b,)E{b,lb,), 
^'/'A],l,X^\bsb^or'x~''[E{b,/bs)E{b,/b,)r\ 
'/'A],\X''{brb,)-'x''E{b,/b,)E{br/b,), 
'/%',{b,bsb,or'x-'' [E{br/b,)E{b,/bj)]-' , 
'/%l{b,brb,)-'x''E{bs/bs)E{b,/bs), 
'lX{b,b,b,b,,)-'x''' [E{b,/b,)E{b,/b,)]-' , 

(B.43) 
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We are able to conclude (to O(x^), 0{x^^), 0{x^^) or 0{x^^)) that the h 
are equal, so that -Bio = b^^- In terms of the functions (5.100) and (5.101), the 
product of the ten Bethe equations gives 



E{x'yb)E{x'yb)E{x^'/b)E{x^yb)E{x^yb)E{x^yb) 



E{x^'^b)E{x^%)E{x^^b)E{x^%)E{x^^)E{x^'^b) 



N 



LdN 



Tsix%)Tiix'b)Ti{x%)Ti{x''b)Tiix''b)J^six'%) 



gs{i/x%)gi{i/x'b)gi{i/x%)gi{i/x''b)gi{i/x''b)gs{i/x'%) 



which can be written as 



-1 



(B.44) 



N 



E{x%)E{x'^b)E{x%)E{x^^b)E{x^%)E{x^%) 
E{xyb)E{xyb)E{xyb)E{x'^^/b)E{x^yb)E{x^yb) 

E{x^yb)E{x^'^/b)E{x^yb)E{x^yb)E{x'^yb)E{x*yb) 
E{x^^b)E{x^'^b)E{x^%)E{x^^b)E{x^%)E{x^%) 



again with elliptic nome x . 



N 



h&N 



(B.45) 
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